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Abstract. We analyze the end-point region of the photon spectrum in semi-inclusive radiative 
Q ■ decays of very heavy quarkonium (ma^ Aqcd)- The S- and P-wave octet shape functions are 

' calculated. When they are included in the analysis of the photon spectrum of the T(15) system the 

agreement with data becomes excellent. 



>- ; The standard Non-Relativistic QCD (NRQCD) factorization [1] (operator product 

expansion) breaks down at the end-point region of the photon spectrum in semi-inclusive 
radiative decays of heavy quarkonium; this is due to the fact that collinear degrees of 
freedom, that are relevant in this kinematic situation, are not included in NRQCD. A 
proper effective field theory treatment of the end-point region of the spectrum requires, 
thus, the combination of NRQCD with the Soft-Collinear Effective Theory (SCET) [2]. 
With this approach factorization formulas has been derived for this process and the 
o • resummation of Sudakov logarithms has been performed [3, 4]. 

If one considers the initial heavy quarkonium state as a Coulombic state, the octet 
shape functions (that appear in the factorization formulas) can be calculated with a 
combination of potential NRQCD (pNRQCD) [5] and SCET. The relevant contributions 
that need to be calculated are depicted in figure 1. From the calculation of the two 
diagrams in that figure we obtain the 5-wave and P-wave octet shape functions in the 
weak coupling regime, the result is the following: 
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FIGURE 1. Color octet contributions. • represents the color octet S-wave current, ▲ represents the 
color octet P-wave current. The notation for the other vertices is ■>:= J^^. Q^p^ Tr [r 6 B] and ©:= 

-^=fi=xTr [r fe E] . The solid line represents the singlet field, the double line represents the octet field and 
the gluon with a line inside represents a collinear gluon. 
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where 



nC&_ _*± = E K_ X i= 1 Ei= _t (5) 
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These shape function are ultraviolet divergent and need to be regularized and renor- 
malized (see [6] for the details of the procedure); some finite parts coming from a linear 
divergence must also be subtracted to have a smooth connection with the leading order 
NRQCD results. 

Since there is good evidence that the T(15) can be understood as a weak coupling 
bound state, we compare the results with the experimental data for the T — > Xy decay 
[7]. The calculation in pNRQCD+SCET is reliable in the region z G [0.7,0.95] (where 
z = 2Ey/M, M being the mass of the heavy quarkonium state), in order to have a reliable 
prediction for the whole spectrum we adapt the interpolation formula used in [3] for the 
inclusion of the color octet contributions, that is 

i dr int _i ^ A W (6) 
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Furthermore the low photon energy region of the spectrum z < 0.4 is dominated by the 
so called fragmentation contributions; these contributions must be added to the previous 
ones to obtain the total decay rate; we will use the formulas presented in section II. B of 
the first reference in [3] and the evaluation of the NRQCD color octet matrix elements 
derived in the appendix of [6] to obtain these fragmentation contributions. 

The comparison with data is depicted in figure 2, the solid curve consists of the 
color octet contribution (with the Sudakov resummation performed in [4] included) and 
the color singlet contribution [3], interpolated according to (6) with the leading order 
NRQCD result (which is reliable away from the end-point), plus the fragmentation 




contributions. No error bars are presented in this plot (neither for the experimental 
data points nor for the theoretical curve) since the aim is just to illustrate that a good 
description of the spectrum is possible (a detailed description of this analysis will be 
presented in [8]). Taking into account that the data points below z = 0.4 have large 
errors we can affirm that the agreement of the theoretical curve with data is excellent. 
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